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ABSTRACT 



We report the detection of pulsed 7-rays for PSRs J0631+1036, J0659+1414, J0742-2822, 
J1420-6048, J1509-5850 and J1718-3825 using the Large Area Telescope (LAT) on board the 
Fermi Gamma-ray Space Telescope (formerly known as GLAST). Although these six pulsars are 
diverse in terms of their spin parameters, they share an important feature: their 7-ray light curves 
are (at least given the current count statistics) single peaked. For two pulsars there are hints for a 
double-peaked structure in the light curves. The shapes of the observed light curves of this group 
of pulsars are discussed in the light of models for which the emission originates from high up in 
the magnetosphere. The observed phases of the 7-ray light curves are, in general, consistent with 
those predicted by high-altitude models, although we speculate that the 7-ray emission of PSR 
J0659+1414, possibly featuring the softest spectrum of all Fermi pulsars coupled with a very low 
efficiency, arises from relatively low down in the magnetosphere. High-quality radio polarization 
data are available showing that all but one have a high degree of linear polarization. This allows 
us to place some constraints on the viewing geometry and aids the comparison of the 7-ray light 
curves with high-energy beam models. 

Subject headings: pulsars: individual (PSRs J0631+1036, J0659+1414, J0742-2822, J1420-6048, JT509- 
5850, J1718-3825) 
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1. Introduction 

The Fermi Gamma-ray Space Telescope (for- 
merly known as GLAST) was successfully launched 
on 2008 June 11. The study and discovery of 7-ray 
pulsars is one of the major goals of this mission. 
Studying pulsars at these high energies is impor- 
tant, because a large fraction of the total available 
spin-down energy loss rate (E = 4ir 2 IPP~ 3 ) is 
emitted in 7-rays. Here / is the moment of iner- 
tia of the star (generally taken to be 10 45 gem 2 ), 
P its spin period and P its spin down rate. By 
studying individual Fermi detections as well as the 
population of 7-ray pulsars as a whole, models for 
the high-energy emission can be constra ined (e.g. 



Harding et all2007t IWatters et al.ll2009h 



The models can be divided into three dif- 
ferent families that place the emitting regions 
at different locations in the pulsar magneto- 
sphere. In the so-ca l led p olar-cap model (e.g. 



Daughertv fc Harding! Il996h the 7-ray photons 



are produced close to the neutron star surface 
(within a few stellar radii) near the magnetic 
axis. At the other extrem e are the outer-gap 



models (e.g. iMorid Il983t ICheng et all Il986t 



Romani fc Yadigaroglul Il995f) . which place the 



emitting region near t he light cylinder. Finally, in 
slot-gap models (e.g. Muslimov fc Harding! 2004 ) 
the particle acceleration occurs in a region bor- 
dering the last open field lines at a large range 
of emission altitudes. The two-pole caustic model 
(jDvks fc RudakH2003l ) is a geometrical realization 
of the slot-gap model. 

It has recently been demonstrated that 7-ray 
pulsars can be discovered via bl ind searches in the 
Fermi data (|Abdo et al.ll2009bl ). Nevertheless the 



detection threshold for pulsed 7-rays is lower when 
accurate positions and spin frequencies (includ- 
ing their unpredictable timing irregularities) are 
already known. Therefore a set of pulsars with 
E > 10 34 ergs -1 is being timed in the radio band, 
allowing Fermi to search f or pulsations with the 
highest possible sensitivity ( Smith et al. 2008). In 
addition these radio observations allow us to de- 
termine the difference in arrival time of the 7-ray 
pulses with respect to the radio pulses, an im- 
portant parameter to distinguish between different 
high-energy models. 

In this paper we report the detection of pulsed 
7-rays for six pulsars that were found by folding 



the Fermi data on the ep hemerides obtained from 
radio observations (e.g. Weltevrede et al.l [2009) . 
These Fermi detections will be included (although 
with more limited c ount statistics) in th e LAT pul- 



sar catalog paper ( Abdo et al. 2009e! ). The 



six 

pulsars of this paper have moderate to large spin- 
down luminosities (E > 10 34 5 ergs -1 , see Table 
[T]) and have (at least given the current count statis- 
tics) 7-ray light curves which are consistent with a 
single peak. Five of the six pulsars have a strong 
degree of linear polarization in the radio band, 
which can be used to constrain the emission ge- 
ometry. Therefore the combination of radio and 
7-ray data for these objects make them valuable 
for tests of the underlying beaming geometry. 

The paper is organised such that we will start 
with describing the Fermi observations in Section 
[2j This is followed by a description of the meth- 
ods used to constrain the emission geometry using 
radio data in Section [3J The results, including 
the 7-ray light curves and spectral parameters ob- 
tained with Fermi, are presented in Section [4] and 
discussed in SectionO Finally the results and con- 
clusions are summarised in Section [6] 

2. Observations 

2.1. Temporal analysis 

The 7-ray data are collected by th e Large 
Area Telescope (LAT: lAtwood et alJl2009h . a pair- 
production telescope on board Fermi. With a 
large effective area (~ 8000 cm 2 above 1 GeV, on 
axis), a broad field of view (~ 2.4 steradian) and 
a high angular resolution, this telescope is far su- 
perior to the Energetic Gamma Ray Experiment 
Telescope (EGRET) on the Compton Gamma Ray 
Observatory (CGRO). The 7-ray events are time 
stamped using a GPS clock on board the satellite. 
These arrival times are transformed to a barycen- 
tric arrival time using the Fermi LAT Science 
ToolB by taking into account the orbit of the 
satellite and the solar system e phemerides (Jet 



Propulsion Laborat ory DE405; IStandish 1998; 



Edwards et al. 2006! ). Tests have shown that the 



resulting prec ision is accurate t o at least a few 
microseconds ( Smith et al. I l2008l) . 

In this paper we use photons which were col- 
lected between 2008 June 30 and 2009 May 22. 



1 http: / / fermi . gsfc. nasa.gov / ssc / data / analysis / scitools / overview, html 
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Table 1: Rotational and derived parameters for six pulsars. 



PSR 


PSR 


P 


E 


T a 


Bs h 


Blc c 


(J2000) 


(B1950) 


(sec) 


(10 35 ergs" 1 ) 


(10 3 V) 


(10 12 G) 


(10 3 G) 


J0631+1036 




0.288 


1.73 


43.6 


5.55 


2.18 


J0659+1414 


B0656+14 


0.385 


0.38 


111 


4.66 


0.77 


J0742-2822 


B0740-28 


0.167 


1.43 


157 


1.69 


3.43 


J1420-6048 




0.068 


104 


13.0 


2.41 


71.3 


J1509-5850 




0.089 


5.15 


154 


9.14 


12.2 


J1718-3825 




0.075 


12.5 


89.5 


1.01 


22.6 



"Spin-down age r c = P/(2P). 

^Magnetic field strength at the surface of the star in Gauss B$ = 3.2 X 10 19 (PP) 1 / 2 . 
c Magnetic field strength at the light cylinder in Gauss B LC = 3.0 X 10 8 (P/P 5 ) 1/2 . 



Only so-called "diff use" class events were used 
( Atwood et al. I l2009h . which are those with the 
highest probability to be caused by 7-rays from 
the source. All 7-rays events with a reconstructed 
zenith angle larger than 105° were ignored to avoid 
the intense 7-ray background caused by cosmic- 
ray interactions in the Earth's atmosphere. 

To produce the 7-ray light curves we used pho- 
tons within an energy-dependent radius 9 < 0.8 x 
^GcV 5 degrees of the pulsar position, requiring a 
radius of at least 0?35, but not larger than the 
6> max of 0?5 or 1?0 as sho wn in Table ffl See the 
LAT pulsar catalog paper (jAbdo et al.ll2009ef ) for 
a detailed discussion of the selection criteria. This 
selection maximizes the signal-to-noise ratio over 
the broad energy range covered by the LAT. In all 
cases the background is estimated from off-pulse 
bins from a 1°— 2° ring around the pulsar using the 
same energy- dependent cut. 

The pulsars discussed in this paper are regu- 
larly observed in the radio band near 1.4 GHz by 
the Fermi Timing Consortium (see Table [3]) . The 
observatories involved are the Parkes 64-m radio 
telescope in Australia, the Lovell 76-m telescope at 
the Jodrell Bank observatory near Manchester in 
England, and the 94-m (equivalent) Nangay radio 
telescope near Orleans, France. The timing pro- 
gram at Parkes is described in IWeltevrede et al 



Table 3: Radio timing parameters. 



( 20091 ) and typ ically involves mont hly observa- 
tions. Nangay (ICognard et al. 2009f ) and Jodrell 
Bank ( Hobbs et al.ll2004) observations are made, 
on average, every 5 to 9 days. 



PSR 


Obs a 


rms 


DM 


(J2000) 




(fis) 


(cm -3 pc) 


J0631+1036 


J/N 


51 


125.36±0.01 


J0659+1414 


P/N 


427 


13.7±0.2 


J0742-2822 


P/J/N 


210 


73.790±0.003 


J 1420-6048 


P 


516 


358.8±0.2 


J1509-5850 


P 


1068 


140.6±0.8 


J1718-3825 


P 


434 


247.88±0.09 



a The radio observatories (P = Parkes, J = Jodrell Bank 
and N = Nancay) involved in the radio timing. 



The times of arrival (TOAs) for PSRs J1420- 
6048, J1509-5850 and J1718-3825 were obtained 
solely from Parkes data. For the other pulsars the 
TOAs obtained from the different telescopes were 
combined before making an ephemeris for the spin 
behavior of the neutron star. The TOAs were 
compared with an initial tim ing solution using 
TEMP02 (jHobbs et alj 120061) producing timing 
residuals which then were fit for the spin-frequency 
and its time derivative, as well as for instrumen- 
tal offsets between data from different observato- 
ries. Most of the pulsars showed strong additional 
deviations in t heir spin behavior, known as tim- 



ing noise (e.g. iHobbs et al.ll2004T ). This timing 



noise is modelled by either adding higher-order 
spin-frequency derivatives or by using the fitwaves 
algorithm within TEMP02. 
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Table 2: Parameters of the analysis of the 7-ray light curves. 



PSR 


n a 
"max 


H b 




<5 C 


FWHM d 


(J2000) 


(°) 










J0631+1036 


1.0 


9 x 10" 


-7 


0.44 ±0.02 ±0.0002 


0.25 


J0659+1414 


1.0 


< 4 x 1(T 


-8 


0.21 ±0.01 ±0.002 


0.20 


J0742-2822 


1.0 


9 x 10~ 


-6 


0.61 ±0.01 ±0.001 


0.10 


J1420-6048 6 


0.5 


< 4 x 1(T 


-8 


0.36 ±0.02 ±0.01 


0.35 


J1509-5850 6 


0.5 


< 4 x 10" 


-8 


0.31 ±0.02 ±0.02 


0.40 


J1718-3825 


0.5 


< 4 x 10" 


-8 


0.42 ±0.01 ±0.006 


0.20 



"Maximum radius from the source to include 7-ray photons for the analysis. 
6 Bin-independent H-test probability that the pulsation would be caused by noise fluctuations. 

c Peak position (phase lag) with respect to the peak of the radio profile (the uncertainties are the statistical and systematic errors 
respectively). 

d Full-width-half-maximum of the peak as a fraction of the pulse period. 

e There is some evidence, although it is statistically not significant given the current count statistics, that the light curves of 
PSRs J1420— 6048 and J1509— 5850 consist of two overlapping components. If this interpretation is correct, then S is 0.26 and 
0.44 for the two peaks of PSR J1420-6048 and 0.18 and 0.39 for PSR J1509-5850. 



The resulting timing model allows an accurate 
assignment of a rotational phase to the 7-ray pho- 
tons (using TEMP02), thereby constructing the 
light curves. The timing parameters used in this 
work will be made available on the servers of the 
Fermi Science Support CenteJl. 

An important parameter necessary to align the 
radio profile with the 7-ray light curve is the 
dispersion-measure (DM), which quantifies the 
frequency-dependent delay of the radio emission 
caused by the interstellar medium. The DM of 
PSR J0631+1036 was measured by comparing the 
TOAs at two widely separated frequencies us- 
ing Jodrell Bank data (0.6 GHz and 1.4 GHz) 
after the templates used were carefully aligned 
so that the components at the two frequencies 
coincide in pulse phase. The DM of the other 
five pulsars were obtained by measuring the de- 
lay across the 256 MHz ban d using Parkes data 
(see IWeltevrede et al.l 120091 and Table OH). The 
systematic error in the alignment of the radio and 
the 7-ray light curves is the combination of the 
DM uncertainty and the rms scatter of the radio 
timing residuals. 

The shapes of the light curves (consisting of all 
photons above 100 MeV) were fitted using Gaus- 
sian functions resulting in peak positions and the 



error bars on the light curves were taken to be 
the square-root of the number of photons in each 
phase bin. The results of the fitting of the light 
curves together with the full-width-half-maxima 
(FWHM) estimated by eye are summarized in Ta- 
ble 03 This tabl e includes the bin-ind ependent H- 
test probability ( de Jager et al. 19891 ) , which mea- 



sures the probability that the observed 7-ray light 
curves are caused by noise fluctuations. 

2.2. Spectral analysis 

The spectral analysis uses the same first 6 
mont hs of Fermi data as the LAT pulsar cata- 
log (|Abdo et al.l l2009eh . The LAT "gtlike" sci- 
enc e tooll performs a m aximum likelihood analy- 
sis (jMattox et al.f l996) to fit phase-averaged spec- 
tra for the six pulsars. Instrument Response Func- 
tions (IRFs) allow proper treatment of the direc- 
tion and energy of each event. We used "Pass 6 
v3" , a post-launch IRF update that addresses inef- 
ficiencies correlated with the trigger rate 3 . Angu- 
lar resolution is poor at low energies: at 100 MeV 
and normal detector incidence, 68% of the photons 
from a point source have reconstructed directions 
within ~ 5° of the true direction, decreasing to 
~ 0.2° at 10 GeV. Therefore, the likelihood analy- 
sis must model not just the pulsar under study, 



2 http: / / fermi . gsfc. nasa.gov / ssc / data / access/lat / ephems / 



3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
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but all neighboring 7-ray sources as well. We 
applied the analysis use d for the 3-month Fermi 
Bright Source List (BSL; lAbdo et al.ll2009df ) with 
an updated model for the Galactic diffuse emis- 
sion to this 6-month data set. Then, as for the 
BSL, we extract events in a circle of radius 10° 
around each pulsar. The likelihood model includes 
all sources up to 17° from each pulsar, with spec- 
tral parameters fixed to the values obtained from 
the BSL analysis for those more than 3° away. 
Spectral parameters for the pulsar, as well as for 
sources within 3°, are left free in the fit. Galac- 
tic diff use emission was m odeled using a GAL- 
PROP (jStrong et al. 2004 ) calculation designated 
54_77Xvarh7S, very similar to that available from 
the Fermi Science Support Center 3 . 

Bri ght 7-ray puls ars like Vela (lAbdo et al 



2009d) or the Crab (|Abdo et alj l2009al ) are ob 



served to have spectra which are well described by 
exponentially cutoff power-law models of the form 



dN 



--— = KE G ^ V exp I — - 

O-E " \ ^cutoff/ 



E 



(1) 



in which the three parameters are the normaliza- 
tion K of the differential flux (in units of ph cm~ 2 
s _1 McV -1 ), the spectral index, T, and the cutoff 
energy, -E C utoff- The energy at which the differen- 
tial flux is defined is arbitrary. We choose 1 GeV 
because it is, for most pulsars, close to the energy 
at which the relative uncertainty on the differential 
flux is minimal. The spectra were fitted separately 
using a power-law and a power-law plus exponen- 
tial cutoff. The difference A of the log(likelihood) 
for the two fits determines the significance <7 C utoff 
for the existence of an energy cutoff, which is de- 
fined to be the test statistic difference \/2A. For 
pulsars with a cr cu toff < 3 the power-law with a 
cutoff did not result in a significantly better fit 
compared to a simple power-law, hence the cutoff 
energy is unconstrained. 

The observed 7-ray energy flux _F b s is the in- 
tegral above 100 MeV of the fitted spectral shape 
times the energy. The luminosity is then 



L 7 = 4irf n F oha D 2 



(2) 



where D is the distance and /o is the flux correc- 
tion f actor which depends on the beaming fraction 
(e.g. IWatters et al.l 120091 ). For outer magneto- 
spheric models fn is thought to be ~ 1, which is 



the value we assume throughout this paper except 
in Table [5] This luminosity can then be compared 
to the spin-down energy loss rate of the pulsar to 
obtain the 7-ray efficiency 



7] = L-y/E. 



(3) 



The results of the spectral analysis are summa- 
rized in Table 3] and the derived luminosities 
and efficiencies can be found in Table [5] Note 
that the distances to all the pulsars except PSR 
J0659+1414 are highly uncertain, as described in 
more detail below. 

Uncertainties on the effective area (< 5% near 
1 GeV, 10% below 0.1 GeV and 20% over 10 GeV) 
and uncertainties in the Galactic diffuse emission 
model dominate the systematic uncertainties on 
the spec tral results as desc ribed in the LAT pulsar 
catalog (jAbdo et al. I l2009ej) . The spectral parame- 
ter uncertainties are ST = (+0.3, —0.1), 5E cuto tf — 
(+20%, -10%), and<5F obs = (+20%, -10%). The 
bias on the integral energy flux is somewhat less 
than that of the integral photon flux, due to the 
weighting by photons in the energy range where 
the effective area uncertainties are smallest. We 
do not sum these uncertainties in quadrature with 
the others, since a change in instrument response 
will tend to shift all spectral parameters similarly. 

3. Deriving emission geometries from ra- 
dio data 

Two important angles used to describe the 
emission geometry of pulsars are the angle a be- 
tween the magnetic axis and the rotation axis and 
the angle C between the line-of-sight and the rota- 
tion axis. A related angle is the impact parameter 
j3 = C — ct, which is the angle between the line-of- 
sight and the magnetic axis at its closest approach. 
These angles can be infer red by applying the rotat- 
ing v ector model (RVM; iRadhakrishnan fc Cookei 
19691 ) to the position angle (PA) of the linear po- 
larization observed in the radio band. This model 
predicts the PA of the linear polarization ip to de- 
pend on the pulse phase <p as 



tan(^> - tp ) 



sin a sin (<j> — tpo) 



sin C cos a — cos £ sin a cos (cj) — <j>o) ' 



(4) 



where ipo an d 4>o ar e the PA and pulse phase cor- 
responding to the intersection of the line-of-sight 
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Table 4: The parameters of the spectral analysis of the 7-ray data 3 



PSR 
(J2000) 


Energy Flux (E > 100 MeV) 
(10 -11 ergcm _2 s _1 ) 


r 


-^cutoff 

(GeV) 


Ccutoff 


J0631+1036 


3.04±0.51 


1.38 ±0.35 


3.6 ± 1.8 


3.2 


J0659+1414 


3.17±0.31 


2.37 ±0.42 


0.7±0.5 


2.6 


J0742-2822 


1.83±0.36 


1.76 ±0.40 


2.0 ± 1.4 


2.0 


J1420-6048 


15.9±2.8 


1.73 ±0.20 


2.7 ± 1.0 


4.6 


J1509-5850 


9.7±1.0 


1.36 ±0.23 


3.5 ± 1.1 


5.1 


J1718-3825 


6.8±1.7 


1.26 ±0.62 


1.3 ±0.6 


4.4 



"The sp ectral analysis is based on the first 6 months of Fermi data as presented in the LAT pulsar catalog paper l|Abdo et al.l 
l2009d) . 

Table 5: The 7-ray luminosity and the efficiency, which are functions of the flux correction factor /n and the 
pulsar distance D. 



PSR 




V 


(J2000) 


(10 35 ergs" 1 ) 





J0659+1414 (0.0032 ± 0.0003)/ n (.D/0.288 kpc) 2 (0.0084 ± 0.0008)/ n (L»/0.288 kpc) 2 

J0742-2822 (0.09 ± 0.02)/ o (L>/2 kpc) 2 (0.06 ± 0.01)/ n (L»/2 kpc) 2 

J1420-6048 (6± l)/n(£>/5.6kpc) 2 (0.06 ± 0.01)/ fi (£)/5.6 kpc) 2 

J1509-5850 (0.75±0.08)/o(L>/2.5kpc) 2 (0.15 ± 0.02)/ n (£>/2.5 kpc) 2 

J1718-3825 (l.l±0.3)/o(£>/3.6kpc) 2 (0.09 ± 0.03)/ n (Z)/3.6 kpc) 2 



with the fiducial plane (the plane containing the 
rotation and magnetic axis) if the emission height 
h cm is small compared to the light cylinder dis- 
tance. In this model the PA-swing is an S-shaped 
curve and its inflection point occurs at (pe- 
lt is found that the degree of linear polarization 
is correlated with E such that virtually all pulsars 
with E > 2x 10 35 er gs" 1 have a linear polarization 
fraction over 50% (Welteyrede fc Johnston! I2 OO8: 
see a lso Qiao et al. 19951: von Hoensbroech et al 



1998; ICrawford et al.ll2001b Uohnston fc Weisbers 
2006). Moreover most of these pulsars have 
smooth PA-swings, making it relatively easy to 
apply the RVM model to the pulsars presented in 
this paper. 

If the emission profile is symmetric around the 
magnetic axis, then the inflection point coincides 
with the middle of the pulse profile. However, co- 
rotation of the emitting region causes the inflec- 
tion point to be delayed with respect to the pulse 



profile. This pulse phase difference A(f> between 
the middle of the profile and the inflection point 
of the PA-s wing can be used to estim ate the emis- 
sion height (jBlaskiewicz et al.l ll991) 



hpA 



Pc 
"87 



A<p = \r lc AcI), 



(5) 



where P is the spin period of the pulsar, c is the 
speed of light and i?LC is the light cylinder ra- 
dius. Because the relative shift of the PA-swing 
wi th respect to the profile is independent of a and 
C ( Dvks et aDl2004h . Equation H can be used to fit 
the PA-swin g even for moderate emission heights 
(|Dvkdl2008l) . Use of these equations ignores ef- 
fects of rotational sweepback of the magnetic field 
lines ( Dvks fc Harding 2004h . pro pagation effect s 
in the pulsar magnetosphere (e.g. Petrova 2006h . 
current-induced distortions of the magnetic field 
and the effects of a finite emiss ion height spread 
and emission height differences (<Dvkdl2008l) . 
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Both the PA-swing and the observed pulse 
width contain information about the geometry of 
the pulsar. Under the assumption that the radio 
beam is symmetric about the magnetic axis, the 
pulse width W (which we take to be the full phase 
range over which we see emission) is related to the 
half opening angle p of the beam via 

cos p = cos a cos Q + sin a sin £ cos (W/2) , (6) 



(|Gil et al.lll98l lLorimer fe Kramerll2005h . There- 
fore, if we know the value of p for a given pulsar, 
Equation[6]can be used as an additional constraint 
on Equation [4] to narrow down the allowed region 
in a — (3 space. The value of p can be estimated 
from the emission height via 



2Pc 



(7) 



assuming a radio beam centred about the mag- 
netic axis whic h is enclosed by the last op en dipole 
field lines (e.g. lLorimer fc Kramerll2005h . In gen- 
eral, although there is some evid ence for cona l 
rings centered on the magneti c axis (|Rankinlll 983) , 
many pulsars are patchy (e.g. iLyne fc M anchester 



1988t lHan fc Manchester! 1200 it iKeith et al.l l2009T 
Weltevrede fc Wrightl I2009T T making both p and 
hpA uncertain. This uncertainty is likely to domi- 
nate the total uncertainty of the radio analysis and 
its impact is discussed for the individual pulsars 
in the next section. Because of poorly understood 
systematics it is impossible to come up with sensi- 
ble error bars on the derived values describing the 
geometry. 

In summary, first we can use the RVM to fit 
the PA-swing (Equation [4]) , resulting in contours 
defining the allowed a — (3 parameter space. Sec- 
ondly, from the offset of the PA-swing with respect 
to the total intensity profile, we can estimate an 
emission height (Equation [5]), which can be trans- 
lated to an opening angle of the radio beam (Equa- 
tion [3 assuming h cm = /ipa). This opening angle 
corresponds to another contour in a — (3 param- 
eter space (Equation [6]), which in the ideal case 
would match with the RVM contours. If not, then 
at least one of the assumptions must be incorrect, 
most likely indicating that the radio beam is asym- 
metric with respect to the magnetic axis. 

Additional constraints on the emission ge- 
ometry can be obtained from measurements of 



the termination shock of the surrounding pul- 
sar wind nebula which can provide a relatively 
m odel-independent est imate of the viewing angle 
£ ( Ng fc Romani 20081 ): unfortunately such mea- 
surements are not yet available for the pulsars 
discussed here. 

4. Results on the individual pulsars 
4.1. PSR J0631+1036 

4-. 1.1. The pulsar and its surroundings 

PSR J0631+1036 was d iscovered as a young 
pulsar bv lZepka et al. ( 19961 ) in a radio search tar- 
geting Einstein IPC X-ray sources. Its DM is very 
high for a pulsar in the Galactic anticenter and 
this is argued to be cau sed by the foregrou nd star- 
forming region 3 Mon ( Zepka et al. 19961 ). In ad- 
dition this pulsar could be interacting with (or be 
embedded in) dark cloud LDN 1605 and therefore 

the distance derived from the DM (3.6±1.3 kpc ac- 

■ ii i 

cording to the lCordes fc Lazioll2002 | model) is pos - 

sibly overestimated. Following IZepka et al. I (|l996l) 
we adopt a distance of 1 kpc consistent with the 
observed X-ray absorption. No pulsar wind neb- 
ula (PWN) h as been found for th is pulsar at radio 
wavelengths (|Gaensler et al. I l2000l) . 



4-1-2. 'y-rays 

PSR J0631+1036 clearly shows 7-ray pulsations 
(see Figure [1} and the light curve features a sin- 
gle broad peak (FWHM 0.25 in rotational phase) 
which lags the radio profile by 0.44 (see Table [5]). 
A power law in combination with an exponential 
cutoff fits the 7-ray spectrum significantly better 
than a single power law, although the cutoff en- 
ergy cannot be accurately determined. 

IZepka et~aT1(ll996h claimed a 7-ray detection by 
EGRET of PSR J0631+1036. The detection was 
very marginal and their light curve does not really 
resemble that seen by Fermi. Their estimated 7- 
ray flux is an order of magnitude larger than that 
obtained from the Fermi data (see Table 0J. 

4-1.3. X-rays 

At the pulsar position a faint ROSAT PSPC 
X-ray source has been foun d, which was to o 
weak to search for pulsations (j Zepka et al. 199G). 
Sinusoidal X-ray pulsations were claimed by 
Torii et al.l (|200ll ) using ASCA data. However, 
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^>1.0GeV 




™ =- 0.1-0.3 GeV 




Nancay 1.4 GHz 




Rotational Phase 



-20 20 

Pulse longitude [deg] 

Fig. 2.— The pulse profile at 1398 MHz observed 
at Nangay (black) of PSR J0631+1036, as well 
as the degree of linear polarization (dashed) and 
circular polarization (dotted). The bottom panel 
shows the PA of the linear polarization (if detected 
above 3a) and an RVM fit. 



Fig. 1. — The Fermi 7-ray light curves of PSR 
J0631+1036 in different energy bands showing two 
full rotational periods with 25 bins per period. 
The photons above 3 GeV are shown in black in 
the second panel from the top. The background 
estimate is shown as a dashed line in the top panel. 
The bottom panel shows the phase-aligned radio 
profile. 

an XMM-Newton observation appears to show 
that the X-ray point source is not associated with 
the pulsar and the p ulsations were not confirmed 



(jKennea et al.ll2002l ). The derived upper limit for 



a X-ray point source is 1.1 x 10 ergs - (0.5-2.0 
keV) for the assumed distance of 1 kpc. This lumi- 
nosity is low compared to the Becker fe Triimpeil 
(Il997h relationship between the X-ray luminosity 
and E, which can be seen as evidence for a larger 
distance to the pulsar. 



4-1-4- Radio 

The pulsar's radio spectrum is relatively flat 
from 1.4 up to at least 6.2 GHz (unpublished 
Parkes datecl ) and its radio profile is highly po- 
larize d (jZepka et al . 1996; IWeltevrede fc Johnston! 
2008). The radio profile consists of four compo- 
nents (see Figure ^ and the outer components 
are strongest at low frequencies. The radio profile 
shows a remarkable deep minimum at the pulse 
phase of the symmetry point. Despite the com- 
plex structure of the profile, it is highly mirror 
symmetric, not only in its shape, but also in the 
spectral indices of the different components. 

We fitted the RVM model (Equation [J) to the 
radio polarization data. The resulting x 2 map is 
shown in gray scale and in the black contours in 
Figure [3] One can see a "banana-shaped" region 



'www. atnf.csiro.au/people/joh414/ppdata 
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100 

a [deg] 



Fig. 3.— The x 2 map of the RVM fit of 
J0631+1036 (in gray scale, a darker color indi- 
cates a lower x 2 )- The lowest reduced x 2 is 1-45 
and the solid contours correspond to reduced x 2 
values that are two, three and four times larger. 
The half opening angles p of the radio beam de- 
rived from the observed pulse width are overplot- 
ted (dotted contours). 



in a — (3 space which contains valid solutions for 
the PA-swing, showing that a is not constrained. 

The remarkable degree of mirror symmetry of 
the radio profile can be seen as evidence that the 
radio beam itself has a high degree of symme- 
try and is centered at the magnetic axis. Un- 
der this assumption the radio emission height 
follows from the offset of the steepest gradient 
of the PA swing and the mirror point of the 
profile (Equation [5|) and is found to be Hpa = 
600 km, which is sligh t ly low er than found by 
Weltevrede fc Johnston ( 20081 ) using data with a 
lower signal-to-noise ratio. Using this emission 
height in Equation [7] we expect a half opening an- 
gle of the beam p of 18° if the emission comes 
from the last open field line. In Figure [3] the con- 
tours of p derived from the observed pulse width 
are overplotted (dotted). One can see that such 
a large value of p suggests that a is close to 90° 
and (3 ~ —4°. If only a fraction of the open field 
lines produces radio emission (i.e. if the beam is 
patchy) then a could be smaller. 



4.2. PSR J0659+1414 (B0656+14) 

4- 2.1. The pulsar and its surroundings 

PSR J0659+1414 was discovered in the sec- 
ond Molonglo pulsar su rvey as a radio pulsar 
( Manchester et al. 19781 ) and this pulsar is the 
slowest rotating and has the lowest E of the pul- 
sars discussed in this paper. The distance of the 
pulsar is well known via parallax measurements 
using very long baseli ne interferometry (288^7 
pc: iBrisken et al.ll2003h . PSR J06 59+1414 is as- 
socia ted with the Monogem ring (jThorsett et al 
2003), a bright diffuse 25° diameter supernova 



remnant easily visible in soft X-ray images of 
the sky. It has a possible pulsa r wind neb- 
ula in optical (IShibanov et al. 200G) and X-rays 
dMarshall fc Schulzll2002l )~ 



4- 2. 2. j-rays 

There was a marginal det ection of pulsed 7-rays 
from this pulsar by EGRET (jRamanamurthv et al. 
19961 ) . The light curve observed by Fermi (see Fig- 
ure 2]) is similar to that seen by EGRET, but the 
7-ray background is much lower in the Fermi data 
because of its superior angular resolution. The 
light curve is single peaked (FWHM 0.20 in rota- 
tional phase) and lags the radio peak by 0.21 in 
phase. 

The pulsar is very weak above 1 GeV, showing 
that its spectrum is extremely soft. Indeed this 
pulsar appears to have the steepest spectrum and 
the lowest cutoff energy of the pulsars in this pa- 
per (see Table 0]). Although the cutoff appears 
to be within the energy range of the LAT detec- 
tor, a power law plus exponential cutoff does not 
describe the data better than a single power law. 
The 7-ray efficiency of PSR J0659+1414 is very 
low (see Table [5]). Because the distance to the pul- 
sar is well determined, the low efficiency cannot be 
caused by an incorrectly estimated luminosity. 

4.2.3. X-rays 

PSR J0659+1414 is one of the br ightest isolated 



neutr on stars in the X-ray sky (jCordova et al 
1989) and it is one of the "Three Musketeers'' 



(the others being Gem inga and PSR B1055-52; 
Becker fc Trumperlll997l ). The X-ray emission is 
a combination of thermal (blackb ody) and non- 
thermal (power law) emission (e.g. iDe Luca et al 
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Rotational Phase 



Fig. 4.— As in Figured] but for PSR J0659+1414. 



20051 ) and is consist ent with a cooling middle -aged 



neutron star (e.g. iBecker fe Trumperl 11997). At 



soft X-rays the pulse fraction is low and the pul- 
sations are sinusoidal (see Figure [5]), typical for 
thermal emission from the surface of a neutron star 
with a non-uniform temperature distribution (i.e. 
hotter polar caps). At higher energies (> 1.5 keV), 
where the non-thermal component dominates, the 
pulsed fraction increases and the profile becomes 
single peaked. 

A relatively aligned geometry is consistent with 
the sinusoidal soft X-ray profile, which suggests 
only one pole is visible from Earth and the low 
amount of modulation in soft X-rays also hint s 
toward an aligned rotator (|De Luca et al. I l2005l) . 
However, there are some complications in mod- 
elling of the soft X-rays, because there is an ap- 





Fermi > 0.1 GeV 



XMM 1.5-7.0 keV 




XMM 0.15-2.00 keV_ 




r-UV HST ~ 



Optical Palomar 
1 1 h- 



Arecibo 1418 MHz 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Rotational phase 

Fig. 5. — The phase-aligned light curves of PSR 
J0659+1414 at multiple wavelengths. The Fermi 
light curve (with its background level indicated 
by the dashe d line) is compared with the XMM 
observations ( De Luca et al. 2005f ). shown in two 
energy ranges encompassing thermal (0.15-2 keV) 
and non-thermal (1.5-7 keV) components (note 
that for the lower panel the vertical axis does not 
start at zero, indicating a large unpulsed fraction 
of the emission) . The near-UV light curve is from 
IShibanov et al. ( 2005f) , who aligned the light curve 
by making use of the similarity with the optical 
light curve. The background - subtra cted optical 
light curve is from Kern et all d2003l) and th e ra- 
dio profile is from Everett fc Weisbergl ([2001). 



parent anti-correlation between the hot and cool 
blackbody component of the thermal part of the 
emission that is not easily understood without 
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Fig. 6.— As in Figure[2j but for PSR J0659+1414 
at 1418 MHz. These A recibo data are taken from 
Weisberg et"aH (|l999h . 



invoking significant multipolc components of the 
magnetic field or magnetospheric reprocessing of 
thermal photons. In addition the best-fitting emit- 
ting radius of ~21 km (using the very accurate ra- 
dio VLBI parallax measurement) is unlikely given 
the expectations for a standard neutron star. 

4.2.4. Optical 

PSR J0659+1414 is seen at optical wavelengths 
( Caraveoetal 1 ll994h . allowing the st udy of the 



(non- thermal) optical pulsations (e.g. iKern et al 
2003). The optical light curve is double peaked 
(see Figure Et and is very simi lar to that seen in 



near- UV ijShibanov et al.ll2005h . The first optical 



peak following the radio pulse is aligned with the 
single 7-ray peak, while the second optical peak 
following the radio peak is aligned with the peak 
seen in (non-thermal) X-rays above 1.5 keV. 

4.2.5. Radio 

The radio profile of PSR J0659+1414 is roughly 
triangular at 1.4 GHz (see Figure with a weak 
shoulder at the trailing edge. The weak shoulder is 



Fig. 7.— As in Figured but for PSR J0659+1414. 
The lowest reduced \ 2 is 0.93. 



associated with weak and broad radio pulses, while 
the radio pulses in the triangle hav e a "spiky" ap- 
pearance ( Weltevrede et al.ll2006b). The strongest 



of thes e radio pulses are argued (jWeltevrede et al 
2006a|) to be similar to transient radio bursts 



seen from the so-called rota t ing r adio transients 
(RRATs; iMcLaughlin et all 120061 ). The radio 
emission is almost completely linearly polarized at 
1.4 GHz, but there is significant depolarization at 
the leading edge. There is also significant circular 
polarization, which is strongest at the trailing edge 
of the profile. Curiously, at 6.2 a nd 8.4 GHz the 



profil e is completely depolarized (jJohnston et al 
120061 and unpublished Parkes data 4 ). 

Figure [7] shows the \ 2 ma P 01 the RVM fit 
for a high signal-to-no ise radio profile taken from 



Weisberg et al. (1999). The results are consistent 



with previous results fe.g. iLvne fc Manc hester 



19881 : lRankinlll993l : lEverett fc Weisberg|l200lh and 



it is immediately clear that, as is often the case, a 
is unconstrained. 

The steepest gradient of the PA- swing lags the 

peak of the radio peak by 14?9±0?7 (jEverett fc Weisberg 
l200ll ), suggesting an emission height of 1200 km 
(Equation [5]) and p — 22° (using Equation [7|) , im- 
plying that a should be ~ 50° (see Figured}. The 
weak shoulder at pulse phases 18°- 30° in Figure 
[5] could indicate that the peak of the profile leads 
the pulse phase corresponding to fiducial plane, 
which would suggest that a < 50°. A relatively 
aligned geometry would be in line with the above 
discussion about the thermal X-rays. The opti- 
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Fig. 8.— As in Figured] but for PSR J0742-2822. 

cal light curve is linearly pol arized and optica l 
PA points can be measured (jKern et al . 2003). 
Unfortunately, while the optical linear polariza- 
tion is potentially useful, the limited data (three 
phase bins) and large errors make it impossible at 
present to use these data to constrain the RVM 
fit. Improved optical polarization measurements, 
however, have the potential to greatly refine our 
geometrical knowledge of this important pulsar. 

4.3. PSR J0742 2822 (B0740 28) 

4- 3.1. The pulsar and its surroundings 

PSR J0742-2822 was discovered as a radio 



pulsar by iFacondi et al.1 (|1973l ). iKoribalski et al 



(1995) determined a kinematic distance between 
2.0 ± 0.6 and 6.9 ± 0.8 kpc, which was higher than 
the distance derived fro m the DM according to the 
Taylor fc Corded 1 19931 mod el (1.9 kpc). This was 
argued by Koribalski et al. ( 19951 ) to be caused by 
an overestimation of the electron density of the 
Gum Nebula in this position. There is a steep 



gradie nt in the electro n dens ity in this direction 
in the ICordes fc Lazid (|2002l ) model, which was 
constructed such that its predicted distance is con- 
sistent with the kinematic distance. We therefore 
adopt a distance of 2 kpc with the note that it 
could possibly be as distant as 7 kpc. 

4-3.2. j -rays 

PSR J0742-2822 has a relatively narrow single 
peak in 7-rays, especially above 1 GeV (see Figure 
[5]). The 7-ray peak lags the radio peak by 0.61 
in phase. In contrast to PSR J0659+1414 this 
pulsar is much weaker at low energies and is not 
detected below 300 McV. Note that if the pulsar 
is at a distance of 7 kpc (the upper limit of the 
kinematic distance), the 7-ray efficiency would be 
~ 70%. This large efficiency may indicate that the 
pulsar is nearer than the distance upper limit. 

4.3.3. Radio 

Like the other pulsars discussed so far, the radio 
emission of this pulsar is highly linearly polarized, 
but there is depolarization at the trailing edge (see 
Figure [9]) . The circular polarization changes sign 
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Fig. 10.— As in Figure H but for PSR J0742- 
2822. The lowest reduced \ 2 is 6.9. 



roughly in the middle of the profile, which is often 
associated wit h emission coming from close to the 
magnetic axis ([Radhakrishnan fe Rankin Il990h . 

The x 2 ma P 01 the RVM fit (Figure [TU]) is much 
better constrained than for the other pulsars. The 
reduced \ 2 of the best fit is not good (6.9), because 
the observed PA-swing deviates slightly from the 
model at the far trailing side of the profile (see 
Figure This deviation is caused by a small 
jump in PA at the pulse phase corresponding to 
the pulse phase where significant depolarization is 
observed. 



Weltevrede fe Johnston! (|2008l ) derived that 
'ipa = 350 km, based on the observed offset 
between the center of the radio pulse profile 
and the location of the steepest gradient of the 
PA-s wing (slightly larger than that d erived by 
e.g. von Hoensbroech fc Xilourid[i~997h . suggest- 
ing p = 18° (Equation [7]). This opening angle 
is inconsistent with the RVM fit (see Figure [TU)) . 
which suggests that the half opening angle of the 
radio beam p < 15° (in line with the value of p 



derived by e.g. iKramer et al. I ll994h . Possibly the 
beam is as ymmetric with respect to the magnetic 
axis (e.g. iLvne fc Manchester 1988f ). Neverthe- 
less, if we believe that most of the open field line 
region is active, this pulsar is unlikely to be aligned 
and we expect /3 ~ —7°. 
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Fig. 11. 
6048. 



As in Figure [TJ but for PSR J1420- 



4.4. PSR J1420 6048 

4-4-1- The pulsar and its surroundings 



PSR J1420-6048 is a 68 ms pulsar (jD'Amico et al 



20011 ) in the northeast wing of the complex of com- 
pact an d extended radio sou rces known as Kook- 
aburra (jRoberts et al.l 120011 ). It is located within 
the 0?32 wide 95% confidence level radius of the 
center of 3EG J1420-6038. The pulsar's large 
spin-down power and distance makes it a plausi- 
ble match for the EGRET source. The situation 
is complicated by the Fermi LAT discovery of pul- 
sations from an X-ray source in the Rabbit pulsar 
wind nebula, PSR J1418-6058 with E = 5 x 10 36 
ergs -1 , only 0?24 away from the r adio pulsar 



and 0?54 from the EGRET source (Abdo et al 



2009b). The period is 110 ms, different from the 
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weak detection of pulsed X-rays with period 108 
ms reported bv iNg et~al ( 2005 ). There are also 
three TeV sources in this region d iscovered by 
HESS (jAharonian et al.l 120061 l2008h . but the as- 
sociation between the TeV sources and the 7-ray 
pulsars is unclear. 

PSR J1420-6048 is the youngest and most en- 
ergetic pulsar of our sample. The pulsar distance, 
derive d from the DM t o be 5. 6 ±0.9 kpc according 
to the Cordes fe Lazio (2002) model is consistent 
with the X-ray absor ption of Ah ~ 2 x 10 22 cm~ 2 
([Roberts et al.ll200ll ) but subject to uncertainties 
because of HII regions and dense clouds in the Ca- 
rina arm. In this paper we apply a more conser- 
vative distance uncertainty of 30% to our adopted 
DM distance estimates (so 5.6 ± 1.9 kpc for this 
pulsar). This is to try to take into account system- 
atic errors in the model for the electron density in 
the Galaxy and this approach i s identical to that i n 
the LAT pulsar catalog paper (|Abdo et al.ll2009eh . 



4-4-2- 'y-rays 

PSR J1420-6048 is the brightest 7-ray pulsar of 
our sample. The light curve of PSR J1420-6048 
(see Figure [XT]) has a broad peak, which probably 
consists of two components lagging the second ra- 
dio peak by 0.26 and 0.44 in phase. This would im- 
ply a peak separation A = 0.18. The second com- 
ponent may follow the common Fermi pulsar pat- 
tern of increasing d ominance at high 7-ray energies 
(jAbdo et al.ll2009eT ) . The small angula r separation 



betwee n this pulsar and J1418-6058 (|Abdo et al 



2009b) considerably increases the background flux 
at the position of the pulsar over a large energy 
interval. The estimation of the background level, 
which follows from a simple measurement of the 
flux in a l°-2° ring around the pulsar, does not 
take into account the flux of J1418-6058, and as 
a consequence the background level is underesti- 
mated. The 7-ray spectrum is fitted significantly 
better by including an exponential cutoff. 

4-4-3- X-rays 

There is a marginal detection of X-ray pulsa- 
tions at the radio pulse period by ASCA (see Fig- 
ure [T2l)_lTOjn_wjU2i^he^C-ray nebula AX J1420.1- 
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6049 ([Roberts et al.l 120011 ). The X-ray and 7-ray 



light curves peak at different phases with respect 
to the radio peak, so it is not clear how the 7-rays 



Fermi > 0.1 GeV 





ASCA (1-10 keV) 
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Fig. 12. — Compilation of phase- aligned light 
curves of PSR J1420-604 8 as seen by Fermi, 
ASCA (jRoberts et alj|200ll) and Parkes. The X- 
ray light curve has an absolute phase error ~ 0.06. 



and X-rays are related. 

4-4-4- Radio 

PSR J1420-6048 is another example of a radio 
puls ar that is nearly completely linearly polarized 
(e.g. Johnston fe Weisberg||2006h and there is also 
some circular polarization (see Figure fl~3|) . The 
radio pulse profile shows a double-peaked struc- 
ture with the trailing component being strongest, 
something that is generally seen for young pul- 
sars with characteristic ag es less than 75 kyr 
(| Johnston fe Weisbergil2006l ). 



RVM mo delling of the PA-swi ng has been car- 
ried out by iRoberts et al. ( 2001 ) who claim a < 
35° and (3 ~ 0?5. It is clear from Figure [Til that 
the best solution can be found at low a values, 
but neither a nor (3 are well constrained. This ge- 
ometry would be consistent with the arc of emis- 
sion around the pulsar seen in the Chandra image 
(|Ng et al.ll2005h . which suggests (if interpreted as 
a torus) a £ which is not particularly small or large. 

IWeltevrede fe Johnstor] (|2008h found that hpA = 
100 km, which is of the same order a s what was 
found by IJohnston fc Weisbergj (|2006l ) (175 km). 
Using the emission height of 100 km in Equation 
[7] leads to an expected p of 15°. This implies (see 



16 




-50 

Pulse longitude [deg] 

Fig. 13. — As in Figure [21 but for Parkes data of 
PSR J1420-6048 at 1369 MHz. 

Figure [T4"|) that the magnetic inclination angle a 
should be relatively small (~ 20°). As always, 
this derivation relies on the assumption that the 
pulsar beam is symmetric and nearly filled. 

PSR J1420-6048 belongs to a group of young 
pulsars with very wide profiles and re latively high 



values of E > 5 x KP 5 ergs (e.g . Manchester 



2005: IWeltevrede fc Johnston! 12008) . They spec- 



ulate that an analogue can be drawn between 
the radio emission of these so-called "energetic 
wide beam pulsars" and their high-energy emis- 
sion. They argued that the sites in the pulsar mag- 
netosphere that produce the radio emission could 
be very similar to those of the high-energy emis- 
sion, leading to the prediction that there should 
be strong similarity between the radio and 7-ray 
light curves. Although the 7-ray light curve may 
indeed be double peaked, the 7-ray light curve is 
significantly offset in phase from the radio profile, 
suggesting a significant difference in the location 
of production of the radio and 7-ray emission. 
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Fig. 14.— As in Figure [3J but for PSR J1420- 
6048. The lowest reduced x 2 is 1-3. 
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Fig. 15.— As in Figure ffl but for PSR J1509- 
5850. 
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Fig. 16. — As in Figure [21 but for Parkes data of 
PSR J1509-5850 at 1369 MHz. 



4.5. PSR J1509 5850 

4- 5.1. The pulsar and its surroundings 

PSR J 1509-5 850 was disco v ered as a 89 ms 
radio pulsar by iKramer et al. ( 2003f h It has a 
pulsar wind ne bula as well as a lon g tail seen 



i n X-ravs (e.g. iKargaltsev et al.l l2008h and radio 
(|Hui fc Beckerll2007D . The distance to this pulsar 
estimated from t he DM is 2.5 ± 0.5 kpc using the 



Cordes fc Lazio d2002h model, but, as discussed 
before, we adopt a more conservative distance un- 
certainty of 30% (2.5 ± 0.8 kpc). 

4-5.2. j-rays 

This pulsar should not be confused with PSR 
B1509-58 (PSR J1513-5908), which was seen in 
soft 7-rays by BATSE OSSE, and COMPTEL 
on CGRO (lUlmer et al.lll993l) . The light curve 
of PSR J1509-5850 (see Figure Q3J) is very broad. 
Similarly to PSR J1420-6048, the light curve may 
be composed of two peaks lagging the radio profile 
by 0.18 and 0.39 in phase respectively (correspond- 
ing to a peak separation A = 0.21). The spectrum 
of this pulsar shows a cutoff at ~ 3 GeV with the 
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Fig. 
3825 



17.— As in Figure ffl but for PSR J1718- 



highest confidence of the pulsars discussed in this 
paper. The derived 7-ray efficiency is very high 
(~ 15%) and is the largest of the pulsars discussed 
in this paper (see Table [5]). However, the distance 
to this pulsar is highly uncertain, therefore the 
luminosity and hence the 7-ray efficiency are not 
well constrained. 

4-5.3. Radio 

This pulsar is the wea kest radio source of our 
sampl e. As noted by IWeltevrede fc Johnston 
(|2008l ). the correlation between E and a high de- 
gree of linear polarization does not hold for this 
pulsar (see Figure [16]). The low degree of lin- 
ear polarization in combination with a low overall 
signal-to-noise ratio prevents us from measuring 
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Pulse longitude [deg] 

Fig. 18. — As in Figure [21 but for Parkes data of 
PSR J1718-3825 at 1369 MHz. 



the PA-swing, and therefore the emission geome- 
try cannot be constrained for this pulsar using the 
RVM model. 

4.6. PSR J1718 3825 

4- 6.1. The pulsar and its surroundings 

PSR J1718-3825 is a 75 m s pulsar discov- 
ered by [Manchester et al~ ( 2001 ) at radio wave- 
lengths. It has a n associated X-ray nebula 



(iHinton et all l2007h an d 



an 



source (jAharonian et al 



associated HESS 
2007J). Its distance de- 



riv ed from the DM is 3.6 ± 0.4 kpc according to 
the lCordes fc Laziol (|2002t ) model, or 3.6 ±1.1 kpc 
by applying a more conservative distance uncer- 
tainty of 30%. 

4-6.2. j-rays 

The light curve of PSR J1718-3825 is single 
peaked (FWHM is 0.20 in phase) and it lags the 
radio profile by 0.42 in phase (see Figure[l7|). Like 
the previous two pulsars discussed its spectrum is 
significantly better described by including a cutoff 
energy. However, the cutoff energy itself cannot 



Fig. 19.— As in Figure H but for PSR J1718- 
3825. The lowest reduced \ 2 is 1-3. 



be determined well. 

4.6.3. Radio 

The radio pulse profile of J1718-3825 has a rel- 
atively complex shape and it is highly linearly po- 
larized (see Figure [T8)) . There is some negative 
circular polarization as well. The 1.4 GHz profile 
is very similar to what is seen at higher frequencies 
(unpublished Parkes data 4 ). 

The PA-swing is S-shaped, allowing an RVM 
fit, although it is not very constraining (see Fig- 
ure [19]). The steepest gradient of the RVM model 
lags the center of the radio pulse by 5°, corre- 
sponding to hpA ~ 80 km (Equation [5J and an 
expected p ~ 13° (Equation [7]) . This would im- 
ply that a ~ 20° (see Figure [T9]) . however, the 
complex shape of the radio profile makes it dif- 
ficult to objectively determine which pulse phase 
corresponds to emission from the fiducial plane. 
If the peak of the profile is a better indicator for 
the magnetic axis than the center, then a is larger 
(~ 50°). 

5. Discussion 

From examination of the light cury es in the 



LAT pulsar catalog (|Abdo et all l2009e[) . certain 



groupings of 7-ray light curves are apparent. One 
group, of wh ich the Vela pulsar is the archetypi- 
cal example (jAbdo et al.ll2009d) . is characterized 
by light curves which consist of two narrow peaks 
separated by ^0.4-0.5 in phase with the first peak 
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offset from the radio profile by ~0. 1-0.2 in phase. 
A second group of light curves also shows dou- 
ble peaks, but the peaks are much closer together 
in phase and can blend together to produce a 
more square-looking light curve. Examples of this 
sort include P SR J1709-4429, the radio-quiet PSR 
J0007+7303 (lAbdo et alJ l2009bh and the radio 
loud PSR J1057-5226 (|Abdo et alJl2009eh . The 
final group of light curves appear to consist of 
a single component, possibly because the peaks 
are completely blended together, are unresolved 
given the present signal-to-noise ratio or simply 
because they are truly single. There is some ev- 
idence that the closer doubles and single compo- 
nent light curves are found preferentially among 
the pulsar with lower E, whereas the Vela-like 
light curves are found f or all spin-down luminosi- 
ties (|Abdo et alJl2009eh . 

The light curves of the six pulsars discussed in 
this paper are all consistent with single peaks. 
Nevertheless, those of PSRs J1420-6048 and 
J 1509-5850 could have the appearance of closely- 
spaced doubles, therefore resembling that of 
PSR J1709-4429. The light curves of the other 
four pulsars do not show, at least given the cur- 
rent count statistics, any hint of a second pulse 
component. PSRs J0631+1036 and J1718-3825 
have relatively broad pulses with an offset to the 
radio pulse of 0.4 in phase. The 7-ray emission 
of PSR J0631+1036 has a tail to earlier phase, 
visible especially above 1 GeV. The light curve 
of PSR J0742-2822 is narrower, and is offset 
from the radio peak by 0.61 in phase. Finally, 
PSR J0659+1414 has a single broad peak at the 
same phase as a corresponding peak in the optical 
and UV and yet lacks the second peak seen at 
those energies. 

We now consider if outer magnetospheric mod- 
els can explain these light curves in conjunction 
with the constraints on the geometry derived from 
the radio profiles by comparing the observed 7- 
ray light curves to those predicted i n the two-pole 



causti c and outer-gap models of IWatters et al 
(2009) 



In the context of the W atters et al. I (|2009h mod- 
els the light curves of PSRs J0631+1036, J1420- 
6048, J1509-5850 and J1718-3825 can all be de- 
scribed in term of closely-spaced double peaks. 
For PSRs J1420-6048 and J1509-5850 there is 
some evidence for a double-peaked nature, while 



for PSRs J0631+1036 and J1718-3825 it is well 
possible that the peaks are unresolved due to a 
lack of signal-to-noise ratio. Double-peaked light 
curves with small separations occur naturally at 
the correct phases relative to the radio emission in 
the outer-gap and two-pole caustic pictures, if one 
skims the edge of the outer magnetosphere cone. 
For narrow gaps (low-efficiency pulsars) this oc- 
curs for a range a ~ 40 — 50° and ( ~ 50 — 65° 
for the outer-gap model and a ~ 45 — 60° and 
C ~ 35 — 55° in the two-pole caustic model. For 
wider gaps (high-efficiency pulsars), the outer-gap 
model has many similar solutions extending to 
a, C ~ 80°, but such pulses are more difficult to 
realize in the two-pole caustic geometry. 

Single peaks that lag the radio pulse by 0.4 - 
0.5 in phase occur in both outer-gap and two-pole 
caustic models for the smallest gap widths, over a 
band spanning a « 55 - 75°, C ~ 20 - 50°. These 
are effectively the second peak of the more typi- 
cal double profile, with the first peak weak or ab- 
sent in this angle range. The preferred geometry 
derived from the radio data of PSR J1420-6048 
(a ~ 20° and j3 ~ 5°) would favour the two-pole 
caustic model (the outer-gap model does not pre- 
dict significant 7-ray pulsations for such a geom- 
etry). However, in this case we would expect the 
7-ray peak to appear at phase 0.1 — 0.2. In any 
case, at present the constraints derived from the 
radio data are too weak to draw firm conclusions. 

If the light curve of PSR J0742-2822 consists 
only of a single narrow peak at phase 0.61 it can- 
not easily be explained by either the outer-gap or 
two-pole caustic model. However, it is tempting to 
associate this peak with the second peak seen in 
the Vela-like light curves as it has the correct phase 
offset with respect to the radio profile. If we do 
this then the first peak must either be missing en- 
tirely or at least be much weaker than the second 
peak; we note that the weak excess at phase 0.12 
in Figure [5] is at the expected phase - longer inte- 
grations should eventually settle this. This picture 
would tie in nicely with the general trend of a rel- 
atively strong second peak at higher energies as 
seen for double-peaked light curves. 

Finally we discuss the curious case of PSR 
J0659+1414. It has a high signal-to- noise ratio 
light curve and is clearly single. The spin param- 
eters of PSR J0659+1414 are not exceptional. Its 
E is relatively low, but that of PSR J1057-5226, 
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which has a harder spectrum, is lower. However, 
the magnetic field strength at the light cylinder is 
very low for J0659+1414. A special geometry may 
be responsible for the unusual properties of this 
pulsar and we consider two possible interpreta- 
tions. In the first, PSR J0659+1414 is more Vela- 
like. In this scenario, the second component (vis- 
ible at lower energies) is entirely missing at 7-ray 
energies. However, although the phase of the visi- 
ble component can then be well explained by both 
outer-gap and two-pole caustic models, it implies 
relatively large values for a and £ which is some- 
what at odds with the constraints from the mod- 
elling of the radio data and the thermal X-rays. 
In the second interpretation PSR J0659+1414 is 
an aligned rotator (a < 40°). In this case in 
the outer-gap picture one would not expect strong 
emission, but the two-pole caustic model does in- 
deed show a single pulse at phase 0.1 — 0.2 later 
than the radio emission. We also note that the 
pulsar has an extremely low efficiency (less than 
~1%) and as its distance is well known and the 
flux correction factor /n is likely to be of order 
unity there is little way to avoid this conclusion. 
Further, its spectrum is amon g the softest of al l 
the pulsars detected by Fermi (jAbdo et al.ll2009ei ) 
and it is virtually undetected at energies > 1 GeV. 
This might also point to a special, aligned geome- 
try where 7-ray emission at somewhat lower alti- 
tude (> 0.1i?Lc) is being observed. 

In summary, therefore, we find that both 
the two-pole caustic model and the outer-gap 
model do a good job in predicting the phase 
offset with respect to the radio emission for 
PSRs J0631+1036, J1420-6048, J1509-5850 and 
J1718-3825 under the assumption that these are 
closely-spaced double-peaked 7-ray emitters. We 
note that the light curve shape, especially of 
PSR J1718-3825, is not well predicted by the 
models. However, it s hould be stress e d tha t the 
model calculations in Watters et al. (2009) are 
idealized and therefore one would indeed only ex- 
pect approximate agreement between the models 
and the data. One could speculate that to obtain 
more realistic physical models one should include 
a larger range of field lines from which 7-ray pho- 
tons are produced. This smooths out the sharp 
peaks in the light curves caused by caustics; these 
sharp peaks are not observed for the pulsars in 
this paper. This also has the effect of making the 



so-called "bridge" of emission between the peaks 
stronger, causing the p eaks in the light cur ve to 



blended together. In IWatters et alT ([2009) the 



two-pole caustic and outer-gap widths increase 
with E, but the emission region is an infinitely 
thin strip on the gap inner edge. For PSR J0742- 
2822, we surmise that there is a "missing" 7-ray 
peak at phase offset 0.1 (and possibly the second 
peak is "missing" for PSR J0659+1414). Without 
speculating about a missing peak, it is hard to 
understand the light curve of PSR J0659+1414 
in the context of outer-gap models. It may be a 
roughly aligned rotator with 7-ray emission from 
lower down in the magnetosphere. 

6. Conclusions 

We report here on the detection of pulsed 7- 
rays by Fermi for PSRs J0631+1036, J0659+1414, 
J0742-2822, J1420-6048, J1509-5850 and J1718- 
3825. These six pulsars are young to middle-aged 
and, except for PSR J1420-6048, have relatively 
small values of E compared to other known 7-ray 
pulsars. In all cases the 7-ray light curves ap- 
pear single peaked (at least with the present count 
statistics), but there is a hint that at least two of 
them have closely-spaced double peaks. As Fermi 
continues its all-sky survey, the quality of the light 
curves will increase, helping to resolve this issue. 

We present high quality radio polarization pro- 
files for these pulsars and discuss their geometries 
in the context of RVM fitting and simple beam 
modelling. Unfortunately, the narrow phase range 
of the radio emission generally leaves a strong de- 
generacy between the fit a and (3 values. This, in 
combination with the limited 7-ray count statis- 
tics makes it difficult to distinguish between single- 
peaked and double-peaked light curves and hence 
to make the comparison with the model predic- 
tions. We show that models where the 7-ray emis- 
sion occurs at relatively high altitudes in the pul- 
sar magnetosphere, such as the outer-gap or two- 
pole caustic models, do a good job in predicting 
the phase of the 7-ray emission relative to the ra- 
dio emission. However, the shape of the 7-ray 
light curve is less well modelled, and additional in- 
puts to the model are needed to explain the strong 
bridge emission in some pulsars and the strong sin- 
gle 7-ray component in PSR J0742-2822 in partic- 
ular. Finally, PSR J0659+1414 warrants further 



21 



attention. It has a peculiar light curve and phase 
offset with the radio profile, its 7-ray efficiency 
is low and its 7-ray spectrum is extremely soft. 
It may be an aligned rotator with 7-ray emission 
arising relatively low down in the magnetosphere. 
Higher signal-to-noise ratio 7-ray light curves in 
combination with possible additional geometrical 
constraints, such as from PWN imaging, will re- 
sult in stronger constraints on the models. 
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